MODELING TOTAL EXPENDITURE ON WARRANTY CLAIMS 
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Abstract. We approximate the distribution of total expenditure of a retail company over warranty 
claims incurred in a fixed period [0, T], say the following quarter. We consider two kinds of warranty 
policies, namely, the non-renewing free replacement warranty policy and the non-renewing pro-rata 
warranty policy. Our approximation holds under modest assumptions on the distribution of the sales 
process of the warranted item and the nature of arrivals of warranty claims. We propose a method 
of using historical data to statistically estimate the parameters of the approximate distribution. 
Our methodology is applied to the warranty claims data from a large car manufacturer for a single 
car model and model year. 



1. Introduction 



Suppose, a retail company sells items each of which is covered by a warranty for a period W. So, 
the company estimates future warranty costs over a fixed period [0, T], say the following quarter, 
based on historical data on sales and warranty claims. Typically, the length of the period for which 
we estimate total warranty cost, i.e. T, if the period we are considering is [0, T], is much smaller 
than the warranty period W. For example, for a car company, usually the warranty period W is 
three years whereas T is a quarter. We assume 2T < W. 

We consider two kinds of warranty policies, namely, the non-renewing free replacement warranty 
policy and the non-renewing pro-rata warranty policy. Under the first policy, the retail company 
agrees to repair or replace the item in case of a failure within the warranty period W . Under the 
second policy, the retail company refunds a fraction of the purchase price if the item fails within 
the warranty period W . The fraction depends on the lifetime of the item and is applicable to 



non-repairable items such as automobile tires; see Blischke and Murthy [1994 page 171, 133] 



The role and importance of warranty costs in the retail industry has increased considerably and a 



considerable amount of research estimates warranty costs; see 


Amato and Anderson 


|1976|, Ja et al. 


2002 




Kalbfleisch et al. 




Kulkarni and Resnick 


20081, 


Majeske 


[2003 . 


Sahin and Polatoglu 


1995 




For the non-renewing free replacement policy case, under highly structured assumptions on 



the sales process and the times of claims, Kulkarni and Resnick |2008| found a closed form expression 



for the Laplace transform of the total warranty cost for a quarter, allowing computation of quantiles. 
Rather than attempting a closed form solution of the Laplace transform, we study approximations 
of the distribution of total warranty cost in a quarter under fairly modest assumptions on the 
distribution of the sales process of the warranted item and the nature of arrivals of warranty claims. 
Depending on the distribution of the cost of individual claims, we approximate the distribution of 
total warranty cost by a normal or a stable distribution. Computation of quantiles by our method 
is relatively straightforward. In the case of companies issuing non-renewing pro-rata warranties, 
we approximate the distribution of the total warranty cost by a normal distribution. 

The advantage in approximating total warranty cost using our asymptotic results is that our 
method does not require strong assumptions on the sales process distribution or on the nature of 
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arrival of claims, and hence is robust against model error. In practice, the times of sales may not 
fit the renewal or Poisson process models; see Section [6j Similar problems are faced when modeling 
times of claims and here also our method based on asymptotic results provides an alternative by 
doing away with the strict assumptions on the distribution of times of claims. 

We discuss methods of estimating the parameters of the normal or stable distribution, which 
approximate the distribution of the total warranty cost in [0, T\. We apply our methods to the 
sales and warranty claims data from a large car manufacturer for a single car model and model 
year. 



1.1. Outline. The following sections are designed as follows. Section \l.2\ reviews some notation. 
In Section [2j we discuss the case of non-renewing free replacement warranty policy. Section [3] 
discusses the case of the non-renewing pro-rata policy. Both kinds of warranty policies use the 
same assumption on the distribution of the sales process. In Section [4j we show that many common 
models for sales processes satisfy our assumptions. In Section [5] we propose a method of estimating 
parameters of the approximate distribution of the total warranty cost in [0, T\. Section [6] applies our 
methods to the sales and warranty claims data from a large car manufacturer for a single car model 
and a single model year. The paper closes with some concluding remarks about the applicability of 
our results and possible future directions. The proofs of the main results are deferred to Section [Sj 

1.2. Notation. The point measure on K C R corresponding to the point x is given by e x , i.e. for 
any Borel set A C K, 



e x {A) 



1 if x £ A, 
otherwise. 



The set of all Radon point measures on K C R is denoted M P {K). Similarly, the set of all non- 
negative Radon measures on K C K is denoted M+(K). 

The set of right continuous functions with left limits from [— W, T] to R is denoted D([— W, T]) 



and the set of continuous functions from [— W, T] to R is denoted C([— W, T]) Billi ngsley] 1999, 
page 80, 121]. Endow D([-W,T]) with the Skorohod topology and C([-W,T]) with the uniform 
topology. 

The set of all one dimensional regularly varying functions with exponent of variation p is written 



RV P [Resnick 



2007, page 24]. Also, we denote conditional expectation of Y conditioned on X as 
E X [Y], i.e. E X \Y\ = E[Y\X}. 

For easy reference, we give a glossary of notation in Section [9} 

2. Non-renewable free replacement warranty policy 



The free replacement policy is the most widely used warranty Blischke and Murthy, 1994, page 



131] and is used for items such as cars, consumer electronics, etc. Under this warranty, the retail 



company repairs or replaces the item in case of a failure within the warranty period W Blischke 



and Murthy, 1994 page 133]. Typically, such policies are non-renewing. 

Sales process: If a warranty claim for an item comes in the period [0,T], the item must be sold 
during the period [— W, T\. As a setting for our approximation procedure, imagine a family of 
models indexed by n. Let S™ be the time of sale of the j-th item in the period [— W, T\. The sales 
process N n (-) is the point process 

N n (t) = J>s»([-W,t]) = |0 : -W < Sf < t}\. 
j 



We further assume that N n (-) is a random element of D([— W, T]) [Billingsley 1999, page 121]. 



WARRANTY CLAIMS 



3 



We define a Gaussian process (N°°(t),t G [— W, T]) having continuous paths, so N°°(-) is also 



a random element of C([— W, T]) |Billingsley 1999, page 80]. Existence of such a process can be 
guaranteed by the Kolmogorov's continuity theorem |0ksendal 2003 page 14]. We assume that 



the sales process N n (-), after suitable scaling and centering, converges in distribution as n — > oo to 
the limiting process N°°(-) in D([-W,T}). 

Times of claims measures: In the n-th model, let Cj^ be the time of the i-th claim for the j-th 
item sold, where we start the clock at the time of sale of the j-th item, so that S™ + C^ f is the 
actual claim time. Assume for all j, the points {C^^i = 1,2,- ••} do not cluster. The times of 
claims measure for the j-th item M™(-) is 

(2.1) M?(A)=J2ec?.( A )> AeB({0,W]), 

i 

where B([Q, W]) is the set of all Borel subsets of [0, W] and M™(-) G M p ([0, W]). Recall, W is the 
warranty period and so only claims in [0, W] will be respected. 

In the n-th model, assume the random measures {Mj l (-), j > 1} are independent and identically 
distributed for all j. Moreover, assume the common distribution of the random measures remain 
the same for all n. We denote the generic random measure describing claim times as M(-), i.e. 

M(-) = M™(-) for all j and all n. 

Claim sizes: We assume that the claim amounts are independent of the times of claims measures 
{Mj'(-) : j > 1} and the sales process N n (-) and claim amounts for different claims are independent 
and identically distributed. 

We consolidate detailed assumptions in the following section. 



2.1. Assumptions. 

(1) Suppose, u(-) is a non-decreasing function in D([— W, T]) which is continuous at the points 
T — W and 0. The family of centered and scaled sales processes in [-W, T] converges weakly 
to a continuous path Gaussian process N°°(-) in D([—W,T]); i.e. 

(2.2) ^f^-^WjV^.). 

Denote the mean function of N°°(t) as 6{t) = £'[A roo (t)] and the covariance function as 
7(*,t) = Cov[N 00 {s),N 00 (t)]. 

(2) For each n, the times of claims measures {Mj 1 (-), j > 1} corresponding to different items sold 
are independent and identically distributed and the distribution of {M"(-) : j > 1} remains 
the same for all n. The random measure M(-) denotes a random element of M p ([0, W]) 
whose distribution is the same as the common distribution of {M"(-) : j > l,n > 1}, i.e. 

M(-) = M™(-) for all j and all n. For each n, the random measures {M™(-) : j > 1} are all 
assumed to be independent of the sales process N n (-). 

(3) The random measure M(-) is a Radon measure with no fixed atoms except possibly at 
and W, i.e. for < x < W, P[M({x}) = 0] = 1. 

(4) We assume M(-) satisfies E[M 2 ([0, W])] < oo. 

(5) For each n, the claim amounts for different claims are independent and identically dis- 
tributed. The common distribution of the claim sizes does not change with n. 

(6) For each n, the claim amounts are independent of the times of claims measures {M™(-) : 
j > 1} and the sales process N n (-). 
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2.2. Asymptotic approximation of total warranty cost distribution. In the n-th model, 
denote the total number of claims for the j-th item sold, that arrived in the fixed period [0, T] by 



(2.3) 



R 7 = E e ^+^ l ([0,T])e C n i ([0, W]), 



and the total number of claims in [0, T] as R n : 

(2.4) R n = J2 e s?+c-MT])e cli mW]) 

{j:-W<S™<T} i 



{j:-W<S?<T} 



r>n 



We require some notation to state the results. Let r : [0, W] — > [0, 1] be a non-negative non- 
increasing function such that r(0) = 1. Recall the random measure M(-) defined in Assumption [2] 
of Section 2.1 and denote its expectation by m(-) = E[M(-)]. Then, r(y)M(dy) is a random Radon 
measure on [0, W] with expectation m(-), such that for all Borel sets A of [0, W], 



(2.5) rh(A) = E 

Now, define for x S [— W, T], 



r(y)M{dy) 



(2.6) 



f [0T _ x] r(y)M(dy), if < x < T, 

= <! I[-x t-x] r (y) M ( d y)> it T — W < x < 0, 

J { J w] r(y)M(dy), if —W < x <T -W. 



Note that 5(-) is a random function, whose interpretation depends on the kind of warranty policy. 
In the free replacement warranty policy, where r = 1, 5{x) gives the number of claims in [0, T] 
for an item sold at time x, i.e. P[5(x) G •] = P[R r { G -I5 1 " = x]. The point-wise expectation and 
variance of <5(-) are given by 



(2.7) 



h(x) = E[6(x)} 



m([Q,T-x]), if < x < T, 
m([-x,T- a;]), if T - W < x < 0, 
m{[-x,W]), if -W<x<T-W, 



and 

(2.8) f 2 {x) = Var[5(x)}. 
Now, we define a function x : D([-W,T}) i-> IRl ' 117 ] by 

(2.9) x(aOC") =x(T-tt) -x((-«)-), x£fl([-lf,T]). 



Recall the Gaussian process A r °°(-) in (2.2). The Gaussian random variable Jj ^ x(iV 0O )( / u)m(dti) 
has expectation /2 and variance ct 2 given by 



[0,W] 



EixiN 00 ) (u)]m(du) 



[0,W] 



E ix (N°°)(u)}r(u)m(du), 



a 



[0,W] J[0,W] 



Cov [x (N°°) (u), X (v)} m{du)m{dv) 



(2.10) 



[0,W] J[0,W] 



Cov [x (N 00 ) (u), X (N°°) (v)} r(u)r(v)m(du)m(dv). 
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In the non-renewing free replacement policy, we choose r(t) = 1. Hence, the measure m defined 
in (2.5) coincides with m(-) = E[M(-)]. Similar simplifications occur in the definitions of 5, f\ and 
/2, as defined in (2.6), (2.7) and (2.8) respectively. The random function 5 when r = 1, is 



(2.11) 



M([0,T-x]), ifO<x<T, 
6(x) = { M([-x,T-x]), ifT -W <x <0, 

Af([-a;,W]), if-IU<x<T-PU, 



and the expectation and variance are given by 
(2.12) h{x) =E[8{x)\ = 



m([Q,T-x}), 
m([— x, T — x]), 
m([— x, W]), 



if < x < T, 

if T - W < x < 0, 

if —W < x < T — W , 



and 
(2.13) 



f2(x) = Var[5(x)]. 



From now on, till the end of Section [2j we use S, fi and f 2 to mean these simplified versions of 
them. We define two constants c\ and c 2 as 



(2.14) 



ci 



■W,T] 



fi{x)v(dx), 



C2 



f 2 {x)v(dx) 



where v(-),fi(-) and /2( - ) are given in (2.2), (2.12) and (2.13) respectively. Also, since we chose 
r(t) = 1 for all < t < W, the parameters p, and a 2 defined in (2.10) takes the simplified forms 



[0,W] 



E[x(N°°)( U )]m(du), 



(2.15) 



a 



[ [ Cov[ X (N c 
J[0,W] J[0,W] 



{u),x (v)] m(du)m(dv) 



Theorem 2.1. Under Assumptions [TJ(4] of Section [2. 1[ the total number of claims R n is asymptot- 
ically normal; i.e. 



J\f(fi,C2 + a 2 ), where J\f(a,b) is the normal distribution with 
mean a and variance b, c\ and C2 are given in (2.14) and jl and a 2 are given in (2.15). 



Let, COST n ([0, T]) be the total warranty cost during [0, T] in the n-th model. Let {X{\ be iid 
with common distribution F representing claim sizes in [0, T\. Denote, SUAij = 5^ = i-^Q f° r an 
j > 1. Then, CO5T n ([0, T]) = ZZi X i = SUM Rn . 

The distribution F of claim sizes is modeled as having a finite or infinite variance. Distributions 



having infinite variance are often assumed to have regularly varying tails [Bingham et al. 1987 



page 344]. When F has infinite variance, we assume F = 1 — F £ RV- a , < a < 2. 

The following theorem allows us to approximate the distribution of CO57^ 1 ([0, T]) based on the 
assumption we make about the claim size distribution F. 



Theorem 2.2. Under Assumptions flTlw of Section 2.1, we approximate the total cost as follows: 



'1) Suppose, the claim size distribution F is such that V 
Then, as n — > oo, 



Jx 2 F(dx) - (JxF(dx)Y 



< oo. 



(2.16) 



COST n ([0,T])-n Cl E 



'nV 



Af(^=p,c 1 + I ^{c2 + a 2 )), 



where JV(-, •), ci, C2, ft and a 2 are the same as in Theorem 2.1 and E = J xF{dx). 
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(2) Suppose, the claim size distribution F is such that F(x) G RV- a , 1 < a < 2. Define, 



b(x) 



l 

l-F 



(2.17) 



(x). Then, as n — > oo, 

COST w ([0,r]) - nciE 



cfZ a (l), 



where c\ is the same as in Theorem 2.1 E = J xF(dx) and Z a (-) is an a-stable Levy motion 



with Z a (l) having characteristic function of the form 



(2.18) 



E [exp(rrZ a (l))] = exp 



1 



itx ax 



-a-l 



d:X 



(3) Suppose, the claim size distribution F is such that F(x) £ RV- a , < a < 1. Define, 



b(x) 



l 

l-F 



(x) and e(x) = f^ x ^ xF(dx). Then, as n 



oo, 



(2.19) 



COST n ([0,T}) -nc^ejn) 
b(n) 



Z Q (ci) + 1{q,=i}Ci logci, 



where ci is the same as in Theorem 2.1 and Z a (-) is an a-stable Levy process with Z a {c\) 
having characteristic function of the form 



(2.20) E [exp(irZ a (ci))] = exp 



ci 



l)ax 



-a-l 



dx + 



(3-1) 



3. Non-renewable pro-rata warranty policy 

The non-renewable pro-rata warranty policy is commonly used for consumer durables such as 
automobile batteries and tires [Blischke and Murthy 1994, page 169]. Under this policy, the 
manufacturer pays a fraction of the cost of the item in case of failure within the warranty period 
W. The fraction depends on the lifetime of the item. So, if an item of cost q, fails after time t from 
the date of purchase, the manufacturer pays the amount q(t), where 

" c b r(t) ift<W, 
otherwise, 

where r : [0, W] — > [0, 1] is a non-negative decreasing function with r(0) = 1. We call the function 
r(-) our rebate function. In many situations, the rebate function is taken to be a linear or quadratic 
function of the lifetime of the item; see Blischke and Murthy |1994 page 172]. 

In this section, since there is no repair or replacement, each item sold can have at most one 
warranty claim. So, the times of claims measure M(-) has the additional property that for any 
Borel measurable set A C [— W, T], M(A) can only assume two values, or 1 and Assumption [I] of 
Section 2.1 is always satisfied. 

3.1. Approximation of the distribution of total cost of warranty claims. As before, let 
C05T n ([O, T]) be the total expenditure on warranty claims during the fixed period [0, T] in the 
n-th model. Let be the lifetime of the j-th item sold. Then, 

COSV l ([0,T}) = Yl c 6Kq:i)^ 1 ([0,W])e 5 ^ 1 ([0,T]). 
{j--S"e[-w,T]} 



Recall the definitions of rh,5, /i,/2 and x given in ( |2.5[ ), (2.6), (2.7), (2.8) and (2.9) respectively. In 
the case of pro-rata warranty policy, the function r(-) in the definition of the random function £(•) 
given in (2.6) is the same as the rebate function r(-) defined in (3.1). Here the random function 
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5(x) is interpreted as the proportion of price spent on warranty claims for an item sold at time 
x. Note that the rebate function r(-) is known and the randomness of <5(-) stems solely from 
the random measure M(-). We denote the mean and variance of the Gaussian random variable 



I-W x{^ co )( u ) r ^ l {du) by fi and a 2 respectively. The forms of jl and a 2 are given in (2.10). We also 
define two constants c\ and C2 as 



(3.2) 



ci 



■W,T] 



fi(x)v(dx), 



C2 



W,T] 



f 2 (x)v(dx) 



where v(-),fi(-) and /2(-) are given in (2.2), (2.7) and (2.8) respectively. 



Theorem 3.1. In case of pro-rata warranty policy, under the Assumptions T][3 of Section 2.1 

COST n ([0,T])-nc bCl K .._ _ 2 , 



Cb\jn 



where q, is the price of each warranted item, c\ and C2 are given in (3.2) and fx and a are given in 



(2.10). 



4. Examples of sales processes 

In the earlier two sections, we have considered the free replacement warranty policy and the 
pro-rata warranty policy. In both cases, the assumptions on the distribution of the times of claims 
measure M(-) as given in Assumptions [3] and [4] of Section 2.1 are modest and a vast class of measures 
qualify. In comparison, the assumption on the sales process N n (-) given in (2.2) is stricter. Here 
we list several sales processes satisfying (12 . 21) . 



Example 4.1. Renewal Processes 

Suppose, N(-) is a renewal process on [0, oo), where the common inter-arrival distribution has mean 
(pi and variance <p2- For the n-th model, define the sales process N n {-) as N n {s) = N{n{s + W)) 
for s £ [— W, T] and define B(-) to be the Brownian motion on [0,oo). Then, from (9.4) of |Resnick 
2007, page 293] or Theorem 14.6 of Billingsley 1999, page 154], we get 



(4.1) 



n -N{n{s + W)) 



s + W 



n 



a 



6f 



B{s + W) 



on D([-W,T}). Define, v{s) 
process is a special case. 



-vi- 



and N°°{s) 



P2 

p\' 2 



B{s + W). The homogeneous Poisson 



Example 4.2. Non-homogeneous Poisson Processes 

Suppose, v : [—W,T] — > [0,oo) is a continuous strictly increasing function and N(-) a homogeneous 
Poisson process on [0, oo) with intensity 1. Now, define the sales process N n (-) as N n (-) = N(ni>(-)) 
and define B(-) to be a Brownian motion on [0, oo). Applying (9.4) of Resnick [2007| page 293] in 
the case of N(-), we get 



(4.2) 



n I —N(ns) 
n 



B(s) 



on D([0,oo)). Define the composition function -0 : D([0,oo)) — > D([—W,T]) by ip(x) = x o u, 



and since tp(-) is continuous Whitt, 1980 Theorem 3.1], using the continuous mapping theorem 



Billingsley 1999 page 21] to (4.2), we get 



n I —N(ns) 
n 



i>(B(s)) 
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on D([— W,T]), which implies 



// | -N(nu(s)) - u(s)j => B(u(s)) 

on D([-W,T]). Define, N°°(-) = B(v{-)) and Assumption [1] of Section [271] holds. 
Example 4.3. Doubly Stochastic Poisson Processes 



Define B(-) and N(-) as in Example 4.2 and let -Do be the subset of non- negative non- 

decreasing functions of D([—W,T]). Assume, there exists a sequence of random elements {A n } of 
Dq independent of N(-) and after centering and scaling the sequence converges to a continuous 
Gaussian process A r 2°(-) in D([— W, T]); i.e. 

(4.3) A"(Q - nu(.) ^ 

= N(A n (-)). Using the fact that iV(-) is 



on D([—W,T]). Now, define the sales process N r 



independent of {A n }, (|42|) and (|43|) yield [Billingsley[ | 1999 [ page 25] 

B(.) 



n(lAT(n-)-(-)) 
-^(A«(.)-nK-)) 



on D{[— W, T]) x -£>([— W, T]), where N^ ^) and B(-) are independent of each other. Further, using 
Billingsley| [1999 page 37] and A n (-)/n =4> u(-) gives 



l- (A»(-) - «/(■)) 




on £>([— W, T]) x L>o x Z?([— W, T]) and from the continuous mapping theorem [Billingsley , 1999 
page 21] and Theorem 3.1 of|Whitt| [1980], we get 



j- (JV(A»(.)) - A n (-)) 
4. (A»(-) -«/(■)) 



JVf (.) 



on W, T]) x Z?([— W,T]). Therefore, applying the addition functional, 

(^(A n (-)) " K")) B(v(-)) + AT-(.) 

on Z)([— W,T]) and the processes S(-) and iV^°(-) are independent of each other. With N°°(- 
B(v(-)) + N£°(-), this model satisfies Assumption [l] of Section 2.1 



Assumption (14. 3|) is modest and Examples 14. 1| or 14.21 satisfy (4.3) 



5. Estimation procedure 

An important estimation question is the choice of n. We interpret n as a measure of the volume 
of sales of the warranted item. So, n should depend on the size of the company and the nature 
of the warranted item. For example, we would expect larger n for an ordinary car model than a 
luxury car model. We assume that for the time period we are considering, say [— W, T], n does not 
change. The non-stationarity of the sales process N n (-) in this period is captured by the functions 
i/(-),9(-) and 7(-, •) given in Assumption [T] of Section 2.1 If we are ambitious enough to predict the 
warranty cost on some time period further in future, say [T, 2T], we will assume that n does not 
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change for the entire time period [— W, 2T]. Thus, we assume n does not change for the entire time 
period we consider. Since we interpret n as a measure of the sales volume, and n does not change 
for the entire time period, we choose total sales in our observed sales data, say total sales in the 
time period [— W, 0], for n. 

We discuss estimation methods for both the non-renewing free replacement warranty policy and 
the non-renewing pro-rata warranty policy. 



0, @ or (§? 
with common 



5.1. Free replacement policy. Which version of Theorem 2.2 should we apply: 
The answer depends on the data of claim sizes. We assumed claim sizes are i.i.d. 
distribution function F. A diagnostic for determining whether data comes from a heavy-tailed 
distribution is the QQ plot Resnick, 2007, page 97]. If F = 1 — F € RV- a for some a > 0, we 
expect the QQ plot to be a straight line with slope ^. If we decide F £ RV- a , we estimate a using 

Chapter 4] 



one of the various estimators of a available in the literature [Resnick, 2007 



Depending 

on the value of our estimate of a, we determine which version of Theorem 2.2|to use. If our analysis 
yields that F ^ RV— a , we verify that F has finite variance and use version (|l ) of Theorem 



2.2 



For versions ([!]), @ and (|3j) of Theorem^ the limit relations in ( |2.16[ ), fl2.17| ) and fl2.19| ) have 



different sets of parameters. We proceed case by case to discuss how we estimate parameters in 
each case. 



5.1.1. Estimation of the parameters in the limit relation in Theorem 2.2, version ([I]). We estimate 
six parameters given in (2.16): c\, C2, /i, ct 2 , E and V. We estimate E by the sample mean and V 



by the sample variance of the claim sizes. 

For the rest of the parameters, we first analyze the sales data and estimate the functions z^(-), #(•) 
and 7(-,-), given in Assumption [T] of Section |2~Tj We assume that we have observed sales for the 
period [— W, 0] and have not observed sales for the period [0, T]. 

One parametric approach for estimating v(-), which is adopted in Section [6j assumes that nu(-) 
follows the Bass model |Bass[ 1969 in the time period where we have observed sales, say [— IF, 0]. 



Since the Bass model describes the pattern of sales from the introduction of an item in the market 
|Bass[ 1969| , this approach gets additional justification when we have sales data of the warranted 
item starting from its introduction in the market. The Bass model for total sales by time t, T(i) 
(adjusted for our clock, since we have sales data for the period [— W, 0]) is given by 



T(t) = 



= n 



1 - exp(-C(t + W)) 
{C/B -l)eM-C(t + W))'' 



1+ , 

where n is the total sales in the time period of observed sales, say [— W, 0]. Hence, using the Bass 
model for nu(-) } we get that v(-) must have the form T{t)/n and to estimate z^(-), we have to 
estimate the parameters B and C. Let v'{t) be the density of at t. We minimize the squared 
error 

o 

l2 



mm 

B,C 



£ [N n (t)-N"(t-l)-nv'(t)]' 



-w+i 



to obtain estimates (B,C). Using this procedure, we fit the Bass model to our observed data on 
sales (say, on the time period [— W, 0]) and then extrapolate on some future time period, say 
[0, T], on which we have no sales data. We denote estimated u(-) as £(•)• Our estimation of u{-) is 
free from any distributional assumption on the sales process N n (-). 

Now, we obtain the residuals {r t = n" 1 / 2 (N n (t) - N n (t - 1) - v(t) + v(t - 1)) : t = -W + 
1, -W + 2, • • • , 0}. These residuals act as surrogates for {N°°{t) - N°°(t - 1) : t = -W + 1, -W + 
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2, • • • ,0} (recall the limit relation in (2.2)). We use standard time-series techniques on {rt ■ t = 
-W + 1, • • • , 0} to get estimates of {6(t),j(t, s):t,s= -W + 1, • • • ,0}. 

We assume that N°°(t) — N°°(t — 1) = TTZt + SCtZ{t), where Z[t) is a stationary Gaussian 
process and SCt is function of t which takes only positive values. Note that this is an additional 
assumption we need for estimation purposes. We have not assumed E[Z(t)] = or Var[Z(i)] = 1. 

We first plot the time plot of {rt ■ t = —W + l,--- ,0}. If the time plot looks stationary, 
we are done and assume TTZt = and SCt = 1- Otherwise, we estimate TTZt and SCt- We do 
moving average smoothing on r t to get TTZt, which estimates the trend. We plot absolute values 
of (rt — TTZt) and fit another moving average estimator to it to get SCt- 

We assume {jt = (r t — TTZt) /SCt ■ t = —W + 1, • • • ,0} act as surrogates for the stationary 
process {Z(t) : t = —W + 1, • • ■ ,0}. We estimate the sample mean I, sample variance s 2 and 
sample autocorrelation function c(-) of {jt}- Hence, {0(t) — 9{t — 1) : t = —W + 1, • • • , 0} is 
estimated as 

6(t) - §{t - 1) = TTZt + ISCt, 

and recover {8(t) : t = -W + 1, • • • , 0}. Similarly, {Cov[N°°{t) - N°°(t - 1), N°°{s) - N°°(s - 1)] : 
t = —W + 1, • • • , 0} is estimated as 

(5.1) Cov[N°°{t) -N°°(t- l),N°°(s) -N°°{s- 1)] = SC t SC s s 2 c(t - s). 



From (5.1 ), {j(-, ■) : t,s = —W + 1, • • ■ ,0} can be computed. 

We also require {9(t),j(t, s) : t G [0, T],s G [— W, T]}. The problem in estimating {0(i),7(i, s) : 
t £ [0,T],s G [— W, T]} is that we do not yet have estimates of {TTZt, SCt, c(s) '■ t G [0, T], s G 
[W,T + W]}. To get estimates of {TTZ t ,S'Ct,c(s) : t G [0,T], s G [W, T + W]}, fit a polynomial to 
both {TTZt ■■ t = -W+l, ■■■ , 0} and {log(SC t ) : t = -W + 1, ■■■ ,0}. We use the fitted polynomial 
values to estimate {TTZt, SCt '■ t G [0, T]}. We also assume c(t) = if t > W, since we only have 
data on sales from [— W, 0]. If we have sales data for a longer period, then it is also possible to 
estimate c(t) for t > W. Then, using estimates of {TTZt,SC t ,c{s) : t G [0,T], s G [W,T + W]} 
we obtain estimates of {9(t),j(t, s) : t G [0, T], s G [— W, T]} following a similar procedure as the 
one used to obtain {(9(i),7(i, s) : t,s = —W + 1, • • • ,0}. Thus, we complete our estimation of 
{§(t),j(t,s):t,s = -W,--- ,T}. 

Now analyze the warranty claims data to get an estimate for the distribution of the times 
of claims measure M(-), given in Assumption [2] of Section |2.1| Recall that the times of claims 



2.1 



measure in the n-th model for the j-th item sold is M"(-). Also, by Assumption |2j of Section 
{Mj'(-), j = 1, 2, • • • ,n} are independent and identically distributed with common distribution as 
that of M(-). For each item j in our sales data, we consider its times of claims measure Mj l (-). If 
an item j has no record of claims, then we assume that Mj 1 = 0. We compute {MJ((x — 1, x]) : x = 
0, 1, • • • , W} with the interpretation that for x = 0, M™((x — l,x]) = M™({0}). From the plot of 

| ( x ' n Sj=i M™((x — 1, x])^j : x = 0, 1, • • • , Wj, we infer a functional form of the mean measure 
m(-) = E[M(-)\ = E[M™ (•)]. Getting a functional form of m(-) is useful be caus e to compute \x and 



a 2 , given in ( 2.15[ ), we have to integrate with respect to m(dx); see Section 6.2 for an example. We 



denote estimated m(-) as rh(-). 

Recall the definition of the parameters c\ and C2 given in (2.14). To estimate c\ and C2, we 
need to estimate first the functions {fi(x) : x G [— W, T]} and {faix) ■ x G [— W, T]}. Actually, 
we estimate {fi(x) : x = -W, -W + 1, • • • , T} and {h(x) : x = -W, -W + !,■■■ ,T}, and get 
estimates c\ = J^ w fi(x)u(dx) and £2 = f^ w f2{x)v(dx) using the trapezoid method of integration. 



We estimate {fi(x) 

(5.2) 

and 



(5.3) 



hix) 



= -w,- 



1 y-m 



WARRANTY CLAIMS 

,T} and {f 2 (x) : x = -W, 



ii 



,T} as 



m([0,T-x]), ifO<x<T, 
m([-x,T-x]), ifT-W<x<0, 
rh([-x,W}), if -W<x<T-W, 



,T-s]) 



Mf([-x,T-x]) 



Ai?([-a?,W]) 



h{x) 



h(x) 2 



if < x < T, 

if T - W < x < 0, 

if - W < x < T — W, 



where n is the total number of items sold and M"(-) is the times of claims measure for the j'-th 
item sold in the n-th model. 

Now, we are left with the estimation of ft and a 2 , given in (2.15). To estimate /2 and a 2 , first we 
must estimate { J B[x(A r °°)(n)] : u G [Q,W]} and {Ccw^iV 00 )^), x(N°°)(v)} :u,v£ [0,W]}, where 
N°°(-) is given in (|]2j) and x(-) is defined in (|]9j). We estimate {E[x(N°°)(u)],u = 0,1,--- ,W} 
and {Cow[x(A roo )(u),x(^ roo )(^)] : v = 0, 1, • • • , W} from the estimates of {§(£}, */(t, s) : t,s = 
-W + l, ■■■ ,T} as 



£[ X (AnW] = #(T-u)-fl(u), 



and 



Co^iV 00 ^)), X (JV°°(i;))] = j{T — u,T — v) + j(-u, -v) - t(T ■ 



-v) -j(T-v,-u), 



where {9(t),j(t, s) :t,s = -W + 1, • • • ,T} are estimates of {9(t),-y(t, s) : t,s = -W + 1, • • • , T} 
obtained while analyzing the sales process. The definitions of the functions #(•) and ■) can be 
found in Assumption [T] of Section 2.1 Now, we integrate by the trapezoid method to obtain the 
estimated mean 



E[x(N°°)(u)]m(du) 



and the estimated variance 
12 



a 



[0,W] J[0,W] 



[Q,W] 



Cov [ X (N°°) (u), X (N°°) (v)} m(du)rh(dv). 



This method of estimation is applied to the sales and claims data of a car manufacturer for a 
specific model and model year in Section [6} 

5.1.2. Estimation of the parameters in the limit relation in Theorem 2.2, version @. We estimate 
the parameters ci, E, a, b(n) and the parameters of the stable distribution of Z a (l), where Z a (l) is 
given in (2.17). Estimate c\ and E in the same manner as described in Section 5.1.1 We estimate 
a by one of its estimators |Resnick , 2007, Chapter 4], say the QQ-estimator. There are two ways 
to estimate b(n): 

(1) Use the (1 — ~)-th quantile of the iid data on claim sizes as b(n); or 

(2) Assume the claim size distribution F is close to pareto and use n l ^ a as an estimate of b(n). 
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We adopt the second method of estimating b(ri) when analyzing data i n Section [6} 



For the stable distribution of Z a (l), we follow the parameterization of Samorodnitsky and Taqqu 
[1994 page 5]. From (2.18), we get that the parameters of the distribution of Z a (l) are Samorod- 
nitsky and Taqqu, 1994, page 171]: 



(5.4) 



A* = 0, 



r(2 



a 



a) ira 

cos( — 

1 v 2 



(3 = 1. 



Obtaining an estimate of a from our estimate of a is a simple numerical procedure. 



5.1.3. Estimation of the parameters in the limit relation in Theorem 2.2, version Estimate 
a, say using the QQ estimator. Depending on whether < a < 1 or a = 1, our estimators of 
parameters will be different, but in both cases, we have to estimate the same set of parameters: 
ci, e(n),b(n) and the parameters of the stable distribution of Z a (a), where Z a {c\) is given in (2.19). 



Estimate c\ using the same procedure discussed in Section |5,1.1| 

When < a < 1, we assume that the claim size distribution F is quite close to Pareto and 



a_ I (l-a)/a 



hence use n l / a as an estimate of b(n) and j 

distribution of Z a (c\), we follow the parameterization of |Samorodnitsky and Taqqu 1994 



l) as an estimate of e(n). For the stable 



page 



5]. From (2.20), wc get that the parameters of the distribution of Z a (c\) arc Samorodnitsky and 
Taqqu[ fl994[^age 170]: 



a = ( ciT(l — a) cos( 



7ra, 



I- a V ' x 2 

Computing estimates of \x and a using our estimates of a and c\ is routine. 

If a = 1, we assume again that the claim size distribution F is quite close to Pareto and hence 
use n as an estimate of b(n) and logn as an estimate of e(n). For the stable distribution of Z a (ci), 
we follow the parameterization of Samorodnitsky and Taqqu| |1994 page 5]. From (2.20), we get 
that the parameters of the distribution of Z a [c\) are |Samorodnitsky and Taqqu 1994 page 166]: 



/i = ci 



sm 



z - zl {z<i}]z dz 



a 



= 1. 



ciir 

/o 

Computing estimates of \i and a using our estimates of a and c\ is then routine. 



5.2. Pro-rata policy case. The estimation method in this case is mostly similar to the one 



described in Section 5.1.1 We need to estimate four parameters given in Theorem 3.1 c\,C2,fi and 
a 2 . 

First, observe that in this case, we do not need any data on claim sizes. Given the times of 
claims measures {M™(-) : j = 1, 2, • • • .}, the claim sizes are determined by the function q(-) given 



in (3.1). 



We analyze the sales process in the same manner as described in Section 5.1. 1[ Thus, we ob- 
tain estimates of the mean and covariance functions of {x(N°°)(u) : u = 0, 1, • • • , W} given by 
{E[x(N°°)(u)] : u = 0, • • • , W} and {Cov[x(N°° (u)), x(N°°(v))] : u,v = 0, • • • , W}. We also esti- 
mate the mean times of claims measure m(-) = E[M(-)] following the same methods described in 



Section 5.1.1 We denote the estimate of m(-) as m(-). 

Now, recall the parameters £l and a 2 given in (2.10) and the rebate function r(-) defined in (3.1). 



Following (2.10), we estimate fl and a as 



[0,W] 



E\x{N°°)(u)]rh{du) 



[0,W] 



EixiN 00 ) {u)]r{u)m{du) 
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and 



[0,W] J[0,W] 



Cov [x (N°°) (u), X (v)} m(du)m{dv) 



a 



Cov [x (N°°) (u), X (A r °°) («)] r(u)r(v)m(du)m(dv). 

'[0,W]J[0,W] 

Now, recall the parameters c\ and C2 given in (3.2). To compute c\ and C2, we first have to 
estimate : x = —W, ■ ■ ■ ,T} and {feix) '■ x = —W, • • • , T}, where the functions and /2O) 

are defined in ([2/7]) and We estimate {fi(x) : s = -W, ■■■ ,T} and {/ 2 (s) : x = -W, • • • , T} 

as 

= { f { _ XiT _ x] r(y)m(dy), 

J[~ x ,w] r (y)™( d y)> 



if < x < T, 

if T - W < x < 0, 

if-VF<x<T-W, 



and 



f2(x) 



1 v-vn 
n Z^j=l 

n 2_fj=l 

n Z^j=l 



V-]Ky)JKT(^) 



h(x 
'fi(x) 
hi?) 



if < x < T, 

if T - W < x < 0. 

if -W < x < T — 



W, 



where r(-) is the rebate function given in (3.1), n is the total number of items sold and M™(-) is 

1, the estimate 



the times of claims measure associated with the j-th item sold. Note that if r 
of {f\{x) '■ x = —W, ■ ■ ■ , T} and {fyix) : x = —W, • • • ,T} is the same as in (5.2) and ( |5.3[ ). Now, 
we integrate by the trapezoid method to estimate c\ = |5w h{ x )v(dx) and 62 = J^ w f2(x)u(dx), 
where £(•) is an estimate of v{-) obtained from the analysis of the sales process. The definition of 
v{ 



is given in (2.2). 



6. Example 

We applied our methods to automobile sales and warranty claims data from a large car manufac- 
turer for a single car model and model year. The company warranted each car sold for three years; 
i.e. W = 1096 days. The period for which we are estimating the cost is taken to be a quarter; i.e. 
T = 91 days. This data is the same as the one used by Kulkarni and Resnick 2008 , but we do not 



assume the sales process or the times of claims measure M(-) is Poisson. 



Which version of Theorem 2.2 should we use? To answer this, we analyze the data on claim 
sizes. 

6.1. Analysis of the claim size distribution. The data consists of the vehicle id which identifies 
the car, the date on which a car comes with some claim, the claim id which is unique for each (car, 
claim) pair and the amount of such a claim. 

From the data, a car on a particular day could come with multiple claims. However, from our 



definition in (2.1) of the times of claims measure MJ l (-) associated with the j-th item sold, MJ l (-) 
is a point measure consisting of random points {Cj^, i = 1, 2, • • • }. So, to be consistent with our 
modeling, for each pair (vehicle id, date), we add the costs of all the claims associated with it, i.e. 
if a car with vehicle id V comes with p claims on a particular date D, which cost X±, X2, • • • , X p 
respectively, then we assume that the car V has arrived on date D with a single claim of size 
X\ + X2 + • • • + X p . Thus, we associate the claim size {X\ + X2 + • • • + X p ) to the (vehicle id, date) 
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Table 1. Summary statistics for claim size data. 



Mean 


Variance 


First quartile 


Median 


Third quartile 


47.53 


18273.14 


7.50 


15.91 


42.79 



density.default(x = actclsize) 



500 1000 1500 2000 2500 

N = 49323 Bandwidth = 2.731 




8 10 



quantiles of exponential 



Figure 1. Density and QQ plot (k 
data points = 49323). 



5000) of the claim size distribution (total 



Our processing of the data on claim sizes differs from that of Kulkarni and Resnick 



pair ( V, D) 
12008) . 

We tabulate the estimated mean, variance and quartiles of the claim size distribution in Table 
[T] Since the estimated third quartile is smaller than the mean, we expect power-like tails of the 



distribution of claim sizes. The density plot of the claim size distribution and the QQ plot [Resnick 
|2007 page 97] are shown in Figure [TJ We use the QQ estimator |Resnick 2007 page 97] to obtain 
an estimate of a = 1.52. This suggests using version ^ of Theorem 2.2 



However, the density plot shown in Figure [T] almost vanishes after the threshold 500, which 
suggests that there are few data points which are relatively very large compared to the rest and 
they are heavily influencing the estimate of a. However, there are 459 data points which are bigger 
than 500. So, on one hand, we cannot discard the claim sizes which are bigger than 500 as outliers, 
while on the other hand, a very small proportion of the data (459/49323 = 0.0093) is influencing 
the summary statistic, the QQ plot and the QQ estimate of a. 

We redo the analysis for all the claim sizes which are less than 500. For this case, the summary 
statistics are tabulated in Tableland the density plot and the QQ plot Resnick, 2007, page 97] 



are shown in Figure [2j Although the data still seems to have a power-like tail, our estimate of a 
using the QQ estimator in this case is a = 2.44, which, to our dismay, suggests using version ([I]) of 
Theorem 12. 2[ 

For comparison, we use both the versions and (|2]) of Theorem [2~2] and compare the quantiles of 
the total warranty cost obtained from the two approximations to check robustness of our asymptotic 
approximation against model error. 



WARRANTY CLAIMS 

Table 2. Summary statistics for claim size data of size less than 500. 
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Mean 


Variance 


First quartile 


Median 


Third quartile 


37.38 


3464.91 


7.41 


15.73 


41.41 



density.default(x = actclsize[actclsize < 500]) 



100 200 300 400 

N = 48864 Bandwidth = 2.636 




quantiles of exponential 



Figure 2. Density and QQ plot (k = 5000) of the distribution of claim sizes which 
are less than 500 (total data points = 48864). 



6.2. Analysis of the distribution of the times of claims measure M(-). The sales data is 
needed to compute the times of claims measures {M™(-) : j = 1,2, ■ ■ ■ , n}. The sales data consists 
of the vehicle id which identifies the car and the date on which it was sold and is a record of 34807 
cars sold over a period of 1116 days. 

To analyze the claims data, we apply the technique in Section 5.1.1 To make any estimation 
about the distribution of the times of claims measure M(-), we obtain the data on {Mj 1 ((z — : 
i = 0, 1, • • ■ ,W}, j = 1,2,- •• , n. Recall, W = 1096 days. For each vehicle id j, note its date 
of sale S™ and the dates on which it comes with a claim, say it comes with claims on p dates 
D\ < Di < ■ ■ ■ < Dp. Now, we compute Cfj = Di — Sj",i = 1,2,- •• ,p. Then, we construct 
the measure M"(-) as Ef=i e Cf (")■ m some cases, we found that Cf^ < (claim honored before 
the car is sold), or C|j > W (claim honored after the warranty period). We handled this as 
follows: if C"j < 0, we make it C™- = and if C-j > W, we make it Cfj = W. Thus, we obtain 
Ml) : = 1 {C^=i} : * = 0, 1, ■ • • , W}, j = 1, 2, • • • ,n. 

times of claims measure m(-) = E[M(-)] in a manner similar to 

The plots of {(i,mi((i-l,i]):=iEi=i^((» * = 
1, i]) : i = 1, • • • , W — 1 = 1095} are shown in Figurepl where 
n = 34807 is the total number of cars in our sales data. Clearly, the plot of {{i,rhi{{i — : i = 
0, 1, • • • , 1096} indicates that the measure m(-) has two atoms at and W = 1096. So, we plot 
{(i, mi((i — 1, i]) : i = 1, • • • , 1095} to infer the structure of the mean times of claims measure m(-) 
in the interval (0, W). The linear appearance of {(i, rhi((i — 1, i]) : i = 1, • • • , 1095} as shown in 



{Mm 



We estimate the expectec 



Kulkarni and Resnick 
1096} and 



0,1,-- - ,W 



2008 
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Figure [3] suggests that for < x < 1096, 

m(dx) = (ax + b)dx. 
By integrating, we get for 1 < i < 1095, 

m((i — = ai + b — -. 

From our fitted line over {(i, rhi((i — : i = I, ■ ■ ■ , 1095} as shown in Figure [3j we obtain the 
estimates 



a = -0.8872 x 10 
We estimate m({0}) and m({W}) by 



b- - = 0.1479 x 10~ 2 . 
2 



mi 



1 n i n 

({0}) := -VM"({0}) = 0.1330, mMW}) := - M? ((W - 1, W]) = 0.0420. 
n z — ' J n z -^ 1 J 



3=1 



i=i 



Thus, we estimate the measure m(-) as m(dx) = (dx + b)dx + m\ ({0 }) l{a ;=o} + w-i({W / })l{x=iy}- 
To estimate the parameters in the limit distribution of Theorem 2.2, version ([!]), we use the 
estimators of /i(-) and /2O) suggested in (5.2) and (5.3). 



6.3. Analysis of the distribution of the sales process. We apply the technique explained in 



Section 5.1.1 We assume that nv(-) follows the Bass model Bass, 1969 for the sales period of 1116 



days. We choose n as the to tal sal es in those 1116 days and so, n 
method discussed in Section 5.1.1 to obtain estimates (B, C) 



34807. We use the least squares 



(4.0149 x 10~ 4 , 1.6738 x 10~ 2 ). The 
time plot of daily count of sales with fitted Bass is given in Figure |4j 

The fit of Bass model is even better for 12-day counts of sale as shown in Figure |4} In case of 
12-day counts, we obtain the least square estimates (B,C) = (4.0279 x 10~ 4 , 1.6740 x 10 -2 ), which 
are not too different from the estimates obtained from daily counts. This gives us confidence in 
our estimates (B, C) = (4.0149 x 10~ 4 , 1.6738 x 10 -2 ) obtained from daily counts and we use these 
estimates for the following estimation procedure. 



Recall from Section |5.1.1 , that we have assumed that {Z(t)} is a stationary Gaussian process 



and the centered and scaled residuals {jt} will act as surrogates of {Z(t)}. We show the time plot 
and the normal QQ plot of {jt} in Figure [5} 




Figure 5. Left: Time plot of {j t }, right: QQ plot of {j t }. 



Now, we follow the procedure described in Section 5.1.1 to estimate {0(t), j(t, s)}. 



6.4. Estimation of quantiles. First, we have to decide the time period for which we want to 
approximate the distribution of total warranty cost. We choose two consecutive periods of length 
91 days starting from the last sales in date, i.e. we choose the next two quarters from the last sales 
date. We denote these two quarters as [0, T] and [T,2T], and accordingly adjust our clock. We 
assume that n remains the same for the entire period [—1116, 2 x 91] (recall that we have sales data 
for a period of 1116 days). 

We compute quantiles of total warranty cost using both the stable and normal approximations 
and compare them. 

For approximation using version (normal) of Theorem 2.2 we follow the method described in 



Section [5 .1.1 to obtain estimates of the six parameters: c\, c%, fx, a 2 , E and V. However, note that 



the parameters ci,C2,jl and a 2 depend on the time-period we are considering, that is the estimates 
will be different for time-periods [0, T] and [T, 2T\. Table [3] gives estimates of these parameters for 
the time periods [0,T] and [T,2T]. For both the time periods [0,T] and [T, 2T], we estimate the 
parameters E = 47.53 and V = 18273.14 using estimates from Table [l] 
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Table 3. Estimated parameters ci,C2,/i and a 2 . 



Time-period 


Cl 


C2 




d 2 


[0,T] 


0.0614 


0.0887 


1.0210 


1.5568 


[T, 2T] 


0.0540 


0.0818 


0.8817 


0.9712 



Table 4. Quantiles for the total cost on warranty claims. 





Time period [0, T] 


Time period [T, 2T] 


p 


p-th quantile 
using version (fl) 
of Theorem 2.2 


p-th. quantile 
using version (|2) 
of Theorem 2.2 


p-th quantile 
using version (fl) 
of Theorem 2.2 


p-th quantile 
using version ((2|) 
of Theorem 2.2 


0.50 


110,694.91 


101,448.27 


97,219.87 


89,224.58 


0.75 


119,449.01 


101,791.20 


104,532.99 


89,539.76 


0.80 


121,618.18 


101,897.93 


106,345.11 


89,637.85 


0.85 


124,146.62 


102,040.29 


108,457.35 


89,768.68 


0.90 


127,327.97 


102,258.94 


111,115.03 


89,969.64 


0.95 


132,043.22 


102,723.28 


115,054.12 


90,396.39 


0.99 


140,888.23 


104,857.40 


122,443.19 


92,357.76 



For approximation using version ([2]) (stable) of Theorem 2.2 we follow the method of esti- 



We obtain estimates of c\ for the time-periods [0, T] and 
We estimate the parameter E = 47.53 using estimates from Table [TJ 



mation described in Section 5.1.2 
[T, 2T] from Table [3} 

We estimate a = 1.52 using the QQ estimator (k 

b(n) = n l l a . Using this estimate of a and (5.4), we estimate of the parameters of the distri- 
bution of Z a (l) (following the parametrization of Samorodnitsky and Taqqu |1994 , as described 
in Section 5.1.2) as fi = 0, a = 1.8688 and (3 = 1. We use J. P. Nolan's software available at 



5000) [Resnick"l |2007[ page 97] and obtain 



|http://academic2.american.edu/^jpnolan/stable/stable.html| j to compute the stable quantiles. 

The quantiles of total warranty cost using both the approximations: version and version ([2]) 
of Theorem 2.2 are listed in Table |4j Note that the quantiles of total warranty cost computed 
using version (|l ) of Theorem |2.2| are bigger than those computed using version ^ of Theorem |2.2[ 
but the difference is not huge as one might have expected since version Q is applicable for the 
heavy tailed data whereas version ([I]) is applicable for the light-tailed data. 

We computed the actual number of claims and the total warranty cost for the periods [0, T] and 
[T, 2T] from our data. Though we cannot test the fit of a distribution from a single observation, 
we do some sanity checks to decide how well the approximations work. 



Start with the actual number of claims. The number of claims in [0, T] is R? Q 



2352 and 



the number of claims in [T,2T] is R^ T2 T] = 1516. Let AFr™(-) be the approximation of the dis- 
tribution function of the total number of claims that arrived in period P using Theorem |2.1| 
We compute AF R n (R p ) for both the time periods P = [0,T] and P = [T,2T]. If AFr»(-) 
were the actual distribution function of R p , then AF R n (R p ) would be uniform on [0,1]. Our 



computed AF R n (R p ) values are AFpn (R 



[o,T]^[0,Tp 



0.5381 and AF R n (R 



[T,2T] ^1T,2TP 



0.0029. The 



probability that a Uniform([0, 1]) random variable take a value more extreme than a number a 
is 2min{P[U < a],P[U > a]}, where U ~ Uniform([0, 1]), which, for the numbers {AFpm (Rp) '■ 
P = [0,T], [T,2T]} are 0.9238 and 0.0058 respectively. These values suggest that {AF R n(?) : P = 
[0, T], [T, 2T}} may be reasonable fits for the distributions of {R n P : P = [0, T], [T, 2T]}. 
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Time period 
P 


Actual cost in 
time-period P 
COST n {P) 


Approximation using 
version ([!]) of Theorem 2.2 
^iFp(C0ST n (P)) 


Approximation using 
version ([2]) of Theorem 2.2 
A 2 F P {COST n {P)) 


[0,21 


148,180.60 


0.9981 


0.9998 


[T, 2T] 


98,992.90 


0.5649 


0.9983 



Now, we compute the actual costs for the time periods [0, T] and [T, 2T] , denoted by COST n ([0, T]) 
and COST n ([T, 2T\) respectively. Let A\F P {-) and ^PpO) be the approximate distribution func- 
tions of the total warranty cost using versions ([I]) and ^ of Theorem 2.2 respectively for the 



period P. The computed values of {A i F P (COST n (P)),i = 1,2, P = [0,T], [T,2T]} are noted in 
Table g For computing {A 2 F P (COST n {P)),P = [0, T], [T, 2T}}, we us ed J. P. Nolan's software 
available at http://academic2.american.edu/~jpnolan/stable/stable.html, If AiFp(-) is the actual 
distribution of COST n {P), then A i F P {COST n (P)) would be uniform on [0, 1]. The probabilities 
that a Uniform([0, 1]) random variable take a value more extreme (as explained before in the pre- 
vious paragraph) than A 1 F P (COST n (P)) is greater than that of A 2 F P {COST n {P)) for both the 
periods [0,T] and [T, 2T], which suggest that for our data on sales and warranty claims of cars, the 
approximation of the distribution of total warranty cost using version ^ of Theorem [23] is better 



than the approximation using version Q of Theorem [23] for both the periods [0, T\ and [T, 2T 



Comparing the quantiles of {AiF P (-),i = 1,2,P = [0,T], [T,2T]} given in Table |] with the actual 
costs {COST n {P) : P = [0,T], [T, 2T}} given in Tableg we arrive at the same conclusion. 

7. Concluding remarks 

We have approximated the distribution of the total warranty claims expenses incurred in a fixed 
period. Our assumptions on the distribution of the sales process N n (-) and the times of claims 
measure M(-) are mild and hence our approximation is applicable in a general context. However, 
we have introduced a lot of independence in our modeling. For example, we have assumed that the 
claim sizes are iid, but in practice this may not true. Similarly, the sales process N n (-) and the 
times of claims measures {M"(-)} or the sales process N n (-) and the claim sizes may be dependent. 
We have ignored such dependences, but allowing for dependence might lead to a more realistic 
modeling and better approximation. 

Our estimation procedure is mostly non-parametric and hence generally applicable. However, 



we have assumed a parametric form for z/(-) using the model proposed by Bass 1969 . Estimating 
v{-) non-parametrically might lead to robustness against model error. 

Another issue is the choice of n in our approximation. We interpret n as a measure of sales 
volume. Though total sales in our observed sales data is a natural candidate for n as we have 
argued, it is not the only candidate. Since n plays an important role in the approximation, the 
choice of n might have a significant impact. 

In Section [6] we have demonstrated the applicability of our method. However, a company 
deciding on reserves to cover warranty cost next quarter should use a more complete and carefully 
collected dataset. 

8. PROOFS 



To prove the asymptotic results, first we state and prove two lemmas. 
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Lemma 8.1. Under the Assumptions [TJ [2] and [4] of Section 2.1 for A 6 
[] E S 3 [exp ( iXn' 1 ' 2 [5{S]) - h(S] 

{j:S*e[-W,T]} 

(8.i) = n 

{j:S"G[-W,T}} 



Q71 



l-A 2 (2n)" 1 (^)-/ 1 (5;)) z + o p (l). 



Proof of Lemma 8A_. Using the fact that for all n, | n?=i a « — EHLi &»l — Yn=l \ a i~bi\, for \ai\,\bi\ < 
1, we get 



n 



{j:5™e[-vy,T]} 



£ E 

{j:S™e[-W,T]} 



exp (an"V2 (s(S?)-fi(S?: 



n 

{i ;S"e[-W,T]} 



E 



{j-.s?e[-w,T\} 



l-\ 2 {2n)- l (5{S])-h{S]))' 
exp (iAn- 1 / 2 (5(5J) - /i(S?))) - 1 + A 2 (2n)- 1 (,5(5;) - Z^S?)) 1 



exp 



(W 1 / 2 (<5(S?) - AOS'?))) - 1 - iXn' 1 ' 2 (5(SJ) - h(SJ)) 



+\ 2 {2ny l {5{S™)-h{S])f 



* E 

{j:S"e[-W,T]} 



E S i 



A 3 



Since for all x, S(x) and f\(x) are bounded by M([0, W]) and m([0,W]) respectively, the above 
quantity is bounded by 



E E 

{j:S™e[-H/,T]} 



A 3 



(M([0,W}) + m([0,W})Y 



A 



A 2 



< A 2 A r "([— W, T]) ^_, 



n 



(M([0,W]) + m([0,^])) ; 



n 



A 



(M([0,W})+m([0,W])Y 



(M([0,TU]) + m([0,VT])) 2 



3.2) 4o. 



The convergence in the last step holds by noting that N n {\—W, T])/n — > u([—W, T]) by Assump- 
tion [T] of Section 2.1 and the quantity within the expectation converges to using the dominated 
convergence theorem. To understand how we use the dominated convergence theorem, first note 
that the quantity inside the expectation is dominated by (ilf([0,W]) + m([0, W])) 2 , which has a 
finite expectation by Assumption [4] of Section |2.1[ On the other hand, the quantity inside the ex- 
pectation is also dominated by q^(M([0,W]) + m({0,W])) 3 , which converges to almost surely. 
Hence, using the dominated convergence theorem, we get as n — > oo, 



E 



_|_6\/n 



(M([0,W]) + m([0,W])Y 



A 



(M([0,W})+m([0,W})Y 



0. 



□ 
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The sales process N n (-) is a non-decreasing process on [— W, T] and hence, induces a measure on 
[— W, T]. In the following, we refer to N n (-) to mean both the sales process in D([— W,T]) and the 
measure it induces. It should be clear from the context what we mean by N n (-). The same rule 
of notation holds for the non-decreasing function v{-) defined in (2.2). Now, we state the second 
lemma. 



Lemma 8.2. Under Assumptions [TJ [3] and [4] of Section 2.1 



J.3) 



(1) The integral of (defined in (2.7)) with respect to the centered and scaled sales process 
converges weakly to a Gaussian random variable, i.e. 

( N n -nu \ 



-W,T] 



n 



[o,W] 



x(N°°){u)m{dv,), 



5.4) 



where v(-) and N°°(-) are given in (2.2), the measure fh(-) defined in (2.5), the function 
x(') is defined in ( |2.9[ ) and the Gaussian random variable Jj w , \ (N°°) (u)rh(du) has mean 

ft and variance a 2 given in (2.10). 
(2) The integral of /20 (defined in (2.8)) with respect to the sales process scaled by n converges 
in probability to a constant, i.e. 
1 



n J[-W,T] 



h(x)N n (dx) 



■W,T] 



f 2 (x)v(dx), 



where v{-) is given in (2.2). 

Proof of Lemma\83$ Step 1: We assume 2T < W. So, T < W - T. First, note that, by using 
Fubini's theorem, we get the following three equations 



(8.5) 
(8.6) 



m([0,T-x]) 



(T-W,0) 



[0,T] 



m([— x, T — x]) 



N n -nv 



n 



(dx) 



[0,T] 



N n -nv 



([0,T -u])m(du), 



N n -nv 



n 



(dx) 



(o,r| 



N n -nv 



n 



([-u,0))m(du) 



J, 



+ 



and 
(8.7) 



-W,T-W] 



m([—x, Wj) 



N n -nv 



n 



(dx) 



(T,W-T) 



[W-T,W) 



[W-T,W] 



N n - nv 



n 

N n -nv 



n 



N n -nv 



n 



([— u, T — u])m(du) 
((T -W,T -u])m(du) 

([-u,T -W])m(du). 



Using (8.5), (8.6), ( |8.7| ) and the definition of x from (2.9), we get, 

'N n -nv 



x 



[0,W] 



n 



(u)m(du). 



Step 2: We use the continuous mapping theorem to prove our result. First, define £ : D([—W,T]) \-t 
by 



i.9) 



ax) 



[o,w] 



X(x)(u)m(du), 



xGD([-W,T]), 
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where x(') is defined in (2.9). 

We show that the continuous functions of D([— W,T]) are continuity points of the function 
£ : x I\o,w] x{x){u)fh{du) . The definition of £ is similar to the well-known convolution functions 
[Feller 1971[ page 143]. Suppose, x n — > x in D([—W,T]) in the Skorohod topology and x is 
continuous. Then, x n — > x uniformly in [— W, T] [Billingsley 1999[ page 124] and 



[0,W] 



X{x n ){u)m{du) 



x(x)(u)rh(du) < sup \x(x n )(u) — x( x )(u)\rh([0,W]) 

[0,W] 0<u<W 

< sup 2\x n (u) - x(u)\m([0,W]) -)• 0, 

-W<u<T 



as n — > oo. So, the discontinuity points of £(•) are contained in D([—W,T])\C([—W,T]). 

Since the limiting process N°°(-) G C([— W, T]) and £(•) is continuous on C([— W, T]), by the 
continuous mapping theorem | Billing sley] 1999, page 21], 



on 



i.e., using 



».9), 



(8.10) 



X 



[0,W] 



iV r ' 



£(iV c 



77 



(u)fh(du) 



[Q,W] 



x(N°°)(u)m(du). 



Hence, by (8.8) and (8.10), we have proved part ([!]) of Lemma 8.2 

@ Using Assumption [1] of Section [2~T| we get N n (-)/n on D([-W,T]). Since Skorohod 

convergence implies vague convergence [Jagers] |1972] , we get 



5-11) 



n 



in M+([— W,T]) (see Section [l~2| ) . By Assumption |3] of Section [2~T] and the definition of 5{-) given 
in (2.6), the random function 5{-) is almost surely continuous at all points except at most at T — W 
and 0. Hence, using definition of /2(') in (2.8), Assumption [4] of Section 2.1 and the dominated 
convergence theorem, we get that f2(x) is discontinuous at most at two points, T — W and 0. So, 
f2(x) can be written as 

' f 2)C+ (x) + d u if -W < X <T-W, 
(8.12) h(x) = \ f2,c+( x )+ d 2, ifT-W<x<0, 

k f 2 , c +{x) + d 3 , if0<x<T, 

where f2 )C +{x) is a continuous non-negative function and di,d2,d^ are three constants. Therefore, 
we get 



-W,T] 



f 2 (x)N n (dx) 



1.13) 



- [ f2. c+ {x)N n (dx) + ^N n ([-W, T - W]) 
n J[-W,T] ' n 

+ —N n ((T - W, 0)) + — iV n ([0,T]). 
n n 



By Assumption [T] of Section 2.1, u{-) is continuous at T — W and 0. Hence, using (8.11) we 
get that the last three terms on the right side of (8.13) converge in probability to d\v{\— W, T — 
W]), d 2 v((T - W,Q)) and d 3 u([0,T]) respectively. Also, by fl8.11[) we get 



W,T] 



f2,c+(*)N n (dx) 



W,T] 



f2, c +(x)v(dx). 



WARRANTY CLAIMS 



2.3 



Hence, using (8.12) we get 



- f f 2 (x)N n {dx) = - I f 2 , c+ (x)N n (dx) + ^N n ([-W,T-W}) 
nJ[-w,T] n J[-W,T] n 



;t] 

+ —N n ((T - W, 0)) + — N n ([0,T}) 
n n 

[ f 2 , c+ (x)v(dx) + ^u([-W,T- W]) 

J[-W,T] n 

+ ^ u{{T -W,0)) + ^u([0,T}) 
n n 



W,T] 



f 2 {x)v(dx). 



□ 



Proof of Theorem 2. 1 , Recall that R% denotes the total number of claims in [0, T] for the j-th item 



sold defined in ( |2.3[ ) and R n denotes the total number of claims in [0, T] defined in (2.4). Using 
Assumption [2] of |2 . 1 [ we get 



i.U) 



R%\SJ = 5(Sf). 



Step 1: From Assumption [2] of Section 2.1, we also get that given the sales process N n (-), the 
random variables {R™ : j > 1} are independent. The characteristic function of the centered and 
scaled R n given the sales process N n (-) is for Ael, 



exp ^Ara~ 1/2 I R n 



■W,T] 



h(x)N n (dx) 



exp | i\n~ 1/2 Yl ( R 7~f^ S j)) 
{j-.s™e{-W,T]} 



= n e 

which, using Lemma |8,1[ can be written as 



exp 



an- 1 / 2 (s(sf)-h(s^ 



H E s ? [l-A 2 (2n)- 1 (5(S™)-/ 1 (^)) 2 ] +Op(l) 
{j:S «e[-w,T}} 

J] [l-A 2 (2n)- 1 /2(^)]+o p (l) 

{j:S?£[-W,T]} 



1.15) 



exp 



exp 



Y, -log[l-A 2 (2n)- 1 / 2 (^) 

{j:SV*G[-W,T]} 



+ Op(l) 



■W,T] 



log(l-A 2 (2n)- 1 / 2 (x)) N n (dx) 



+ o P (l) 
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Step 2 : Now, since | — log(l — x) — x\ < 2\x\ 2 if < i, we get that for large enough n, 



(8.16) 



A 2 



2n 



< 



2l 

2n 2 ' 



A 



fi(x)<^[E[M 2 ([0,W})]]- 



C 



where C = ^[E[M 2 ([0,W])]]' z . Hence, from Lemma |8.2| part ^ and fl8.16[ ) it follows that 



/- 



r 



log 1 



2n 



2n 



/ 2 (x) A/ n (cix) 



/. / 2 

/5w ^f2{x)v{dx) 



and so, 



(8.17) 



cxp 



log 1 



2n 



h{x) N n (dx) 



exp 



exp 



T \2 



h{x)v(dx) 



A 2 c 2 



Therefore, using (|8.15|), (|8.17|) and the dominated convergence theorem, we get for A € 

A 2 c 2 



(8.18) 



E 



cxp 



iXn- 1 ' 2 R n 



■W,T] 



A(x)N n (dx) 



exp 



Step 3: Now, we prove the joint convergence of 

X n A _ J_ ( R n -^ 



-W,T] 



fi{x)N n (dx) 



$ { _ WTl h{x)N n (dx)-nc ) ■ 

Consider the joint characteristic function E [exp (i\X n + icj)Y n )] for (A, 4>) £ M 2 and note that 

A 2 c 2 



exp (i\X n + i(j>Y n ) - exp 
A 2 c 2 



cxp 



+ i(j)Y n 



E [exp (i\X n + i<j)Y n )\ = E 
(8.19) +E 

First, we deal with the first term on the right side of (|8.19|). Note that 
E 



+ i(j)Y n 



exp (i\X n + i<f>Y n ) - exp 



A 2 c 2 



+ i(j)Y r 



(8.20) 



= E 
= E 
< E 



exp (i\X n + i<pY n ) - exp 



A 2 c 2 



+ i(pY n 



exp(icf)Y n ) 



E Nn ^ [exp (i\X n )} - exp 



E Nn ^ [exp (i\X n )] - exp 



A 2 c 2 



A 2 c 2 



0, 



where the last convergence follows from (8.18). For the second term on the right side of (8.19), 
observe 



Y h 



-W,T] 



h(x)N n {dx) -nc fy/n 



-W,T] 



N n -nv 



ii 



(dx) 
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Using Lemma 8.2, part Q, we get 

A 2 c 2 



3.21 



lim E 

n— >oo 



exp 



+ i(l)Y r 



exp 



A 2 



C-2 



b 2 a 2 



where the parameters fi and a 2 are given in (2.15). Therefore, using (8.19), (8.20) and (8.21), we 
get 



lim E [exp (i\X n + i<pY n ) 

n— >oo 

So, the joint convergence holds, i.e 



lim E 

n— >oo 



exp 



yn 



A 2 c 2 



x 

Y 



exp 



A 2 c 2 



b 2 a 2 



on M 2 , where X and V are independent, X ~ A/"(0, c 2 ) and Y ~ M(fl,a 2 ). Hence, 

n 

which gives us the required result. 



ci 



□ 



Proof of Theorem 2.2. Let Xi denote the size of the i-th claim which arrived during the time interval 
[0,T]. Denote, SUMj = £f =1 X i for a11 3 > 1 - Then, COST n ([0, T]) = Y%=i X { = SUMi 



(1) By Assumption [5] and Donsker's theorem Billingsley 1999, page 146], we know 

SUM[ n .] - [n-]E 



fnV 



on D([0, 00)), where W(-) is a standard Brownian motion. Also, from Theorem 2.1 we know that 
y/n — ci) converges in distribution to a random variable V, where Y ~ AA(/2, c 2 + <t 2 ). These 
two facts, coupled with Assumption [6] of Section 2.1 gives that 



SUM[ n .]-[n-]E 
n (^T- c l) 



Y 



on D([0,oo)) x M, where W(-) and y are independent of each other. Hence, from Billingsley 1999 
page 37], we get 

SUM< n .]-[n-]E 



(8.22) 



VnV 
R n 



\ n 




on -D([0, 00)) xKxl, where each three components on the right of (8.22) are independent of each 
other. Applying Theorem 3 of Durrett and Resnick [1977] to (8.22) yields 

(8.23) 



SUM 



[R n 



-[R n -]E 



n 



nV 



W( C1 (-)) 
Y 



on D([0, 00)) x M. Since the first component of the limit in ( |8.23 ) is a continuous process, we have 



COST n {[Q,T])-R n E 
- 



W( Cl ) 
Y 
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|Durrett and Resnick 1977, Section 4]. So, 



COST n ([0,T}) -nciE COST n ([0,T]) - R n E E 



'nV 



'nV 



+ 



V V n 



on M., which completes the proof of part ([I]) of Theorem 2.2 

(2) By Assumption[5]and a minor variant of the central limit theorem of heavy-tailed distributions 
|Resnick 2007, page 218], we know 

SUM [M - [n-]E 

a{ ' } 



on D([0,oo)), where Z a {-) is an a-stable Levy motion whose characteristic function satisfies (2.18). 
Now, applying Theorem 3 of Durrett and Resnick 1977 in a similar way as in the proof of (2.16), 
we get 



(8.24) 



5WM [fln .] - [R n -]E 
b(n) 



on D([0,oo)). Since the limit process in ( |8.24[ ) is continuous in probability at time 1, we have 

COST n ([0,T}) - R n E 



(8.25) 
on M. Observe, 
(8.26) 



b(n) 



Z a (ci) =cfZ a (l) 



COST n ([0,T]) - n Cl E COST n ([0,T\) - R n E R n E - nc x E 



b(n) 



b{n) 



+ 



b(n) 



In (8.25), we have already found the weak limit of the first term on the right side of (8.26). 
Now, we deal with the second term on the right side of (8.26). Notice, b(-) G RVi, and hence 

a 

x l / 2 (b{x))~ l € RVi_i. Since a < 2, ^ — ^ < 0. Hence, hnin^oo = 0. Therefore, using this 
fact, along with Theorem 2.1 and Slutsky's theorem, we get 



R n E - n Cl E 



nE 



R n 
n 



< . i 4 o, 



b{n) b(n 

which together with (8.25) and (8.26) completes the proof of part ^ of Theorem |2.2 

(3) By Assumption |5| and the central limit theorem of heavy-tailed distributions |Resnick, 2007 
page 218], we know 

SUM[ n .] - [n-]e(n) 



b{n) 



on D([0, oo)), where Z a (-) is an a-stable Levy process and characteristic function of Z a [c\) satisfies 
( 2.20p . Now, applying Theorem 3 of Durrett and Resnick [1977 in a similar way as in the proof of 
(2.ie|, we get 



5.27) 



SUM 



[R n -]e(R n ) 



b(n) 



Z«(ci(-)) 



on D([0,oo)) and since the limit process in (8.27) is continuous in probability at time 1, we get 
(8.28) COST^-^em ^ ^ 



b(n) 
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on R. Now, notice that 

COST n ([0,T}) -ncfein) 



(8.29) 



COST n ([0,T}) -R n e(R n ) R n e{R n ) - nc? e(n) 



b{n) 



b(n) 



+ 



b{n) 



We have already observed in (8.28) the weak convergence of the first term on the right side of 
(8.29). We now turn to the second term. To show the convergence in probability of the second 
term, we consider two separate cases, namely a = 1 and a < 1. 

a < 1 case: Notice e(x) = I(b(x)) where I(x) = f* tF(dt) £ RVt- a [Resnick[ [2007[ page 36, 
Ex. 2.5]. Also, b(x) £ RVi , and lim. 

page 32], e 



b(x) = oo. Therefore, from Proposition 2.6(iv) of 



2007 



£ RVi_ v So, xe(x) £ RVi. Also, from Theorem 



converges in probability to c\. Hence, |Resnick 2007 page 36] 
(8.30) 



2.1 



Resnick 
we know that R n /n 



R n e{R n ) P 



neinj 



Using (8.30) and the fact that the sequence { i converges to and hence is bounded, we get 
that 



R n e(R n ) -ncfe(n) _ ne(n) / R n e(R 



b{n) 



b{n) 



ne(n) 



i 



0. 



which, together with (8.28) and 

a = 1 case: In this case, we may write the second term of the right side of (|8.29|) as 



3.29) completes the proof for the case a < 1. 
y write the s 

R n e{R n ) - ncie(n) R n e(R n 



b{n) 



3.31) 



b{n) 

R n e(R n ) - e(n) 



R n e(n) R n e(n) — ncie(n) 



b(n) 



n 



b(n)/n 



+ 



neln 



b(n) 



R n 
n 



ci 



First we deal with the first term on the right side of (8.31). Note that F(y) is (-l)-varying 
and so, using a theorem of de Haan 1976] (Proposition 0.11 of Resnick 1987 page 30]), we get 
that Jq F(y)dy is II- varying with auxiliary function xF{x). Hence, using the fact that I{x) = 
Jq F{y)dy — xF(x), it easily follows from the definition of n- varying functions 



Resnick 



1987 



page 



27], that I(x) is also n- varying with auxiliary function xF(x ). Since, b(x) £ RV\, e(x) = I(b(x)) is 
also n-varying with auxiliary function b{x)F{b{x)) = b{x)/x de Haan and Resnick, 1979 , Resnick 



2007[ page 38], 

Using local uniform convergence of n-varying functions on sets away from [Bingham et al. 



1987 



page 139, Theorem 3.1.16], and the fact R n /n — > c\ > 0, we get 
e{R n ) - e(n) _ e (n • f^) - e(n) 



b(n)/n 



Hence, 
(8.32) 



R n e(R r ' 



b(n)/n 

R n e(n) R n e(R n ) - e(n) P 



log Cl. 



b{n) 



n 



b(n) /n 



Cl log Cl. 



Now, we turn to the second te rm in (8.31). Note that b{n) £ RV\ and since e(n) is n-varying, 

ne(ra) 



e(n) £ RVq |Binghamet al 

I 



1987 



page 128]. Therefore, 



b(n) 



£ RV_ i , and so, lim r , 

2 



/ne(n) 



0. 
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2.1 



Also, from Theorem 
Slutsky's theorem, we get 

(8.33) 
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we know that yfn — c±) is asymptotically normal. Therefore, using 



R n e{n) — nc\e[n) 



b(n) 



ne{n) _ / R n \ P 

- r — y/n[ ci I ->-0. 

o{n) \ n 



Finally, using ( |8.3lD , fl8.32| ) and ( |8.33D we get 

R n e(R n ) - recie(n) p 
6(n) * 



Cl log Cl, 



which, together with (8.28) and (8.29) completes the proof for the case a = 1. 



□ 



Proof of Theorem 3.1 Denote by CO5T™([0, T]) the warranty claims expenditures which are in- 
curred in [0, T] for the j-th item, sold in the time interval [— W, T]. In notation, 

COST]([0,T}) = Cfe r(C™ 1 )^ 1 ([0,^])e S7+ ^ 1 ([0,r]). 
So, COST n ([0,T]) = j: {r .s^[-w,T ]} ' 



COST^{[0,T]). Using Assumption [2| of Section 



2.1 



we get 



( Cfe )- 1 CO577([0,T])|S , ™ = ^™ 



Now, we consider the characteristic function of 

(^c b y l cosT n (%T})- 

given the sales process N n (-). Note that for A G 



-W,T] 



h{x)N n {dx) /^n 



E Nn (-) 



exp ( i\n- 1 ' 2 \{c l )' l COST n {[Q,T\) - [ /i(x)iV n (dx 
V V J[-W,T\ 



exp i 



An" 1 ' 2 



c,)- 1 CO577([0,T])-/ 1 (5" 



n 



exp 



{i:5^e[-W,T]} 



Un- 1 / 2 (5(57)-/i(^))) 



Now, we want to use the same steps as in the proof of Theorem 2.1 To do so, we must be able 



to use Lemma 8.1| and Lemma |8.2[ which need Assumptions [T]j4] of Section |2.1[ By the hypothesis 
of Theorem 3.1 



we have Assumptions [I]{3] of Section 2.1 Since M([0, W]) < 1, we know that 
Assumption 4] of Section 2.1 is also satisfied. Hence, we can apply the results of Lemma 8.1 and 



Lemma 8.2 here, too. Now, the proof follows exactly similar steps as in the proof of Theorem 2.1 



So, the rest of the proof is omitted. 



□ 



9. Glossary of notation 



S™ = time of sale of the j'-th item in the n-th model, 
N n (t) = total number of sales in [— W, t] in the n-th model = eg»([— W,t]), 

3 

u(t) = the centering function of ^N n (t), 
A r °°(-) = the Gaussian process limit of centered and scaled -iV"(-), 



M(-) 

m(-) 
mi(-) 
m(-) 



m(-) 
5(z) 

/ 2 (x) 

Cl 
C2 

E 
V 
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= the time of the i-th claim for the j-th item sold, where we start our clock at the time 
of sale 

= e c™ i ( - ) = the times of claims measure for the j-th item sold, 

i 

- the generic random measure representing times of claims for an item sold, i.e. 

M (■) = M™(-) for all j and all n. 
--E[M(-)\, 



n 



estimate of m(-), maybe using a functional form, 

total number of claims in [0, T] for the j-th item sold = > es n +c n - ([0> T])ec n . ([0, W]), 



R n = 



r(t) 



total number of claims in [0,T] = ^ i?", 

{i:S"e[-W,T]} 

' 1, 



if we follow the non-renewing 
free replacement policy, 



the fraction of the cost refunded if the claim if we follow the non-renewing 
comes after time t from the date of sale, pro-rata policy, 



E 



[ r(y)M(dy) 
■/(•) 



' J [0T _ x] r{y)M(dy), if < x < T, 
< j hxT _ x] r(y)M(dy), if T - W < x < 0, 
^ J { _^ w] r(y)M(dy), if —W < x < T — W, 

{m([0,T-x]), ifO<x<T, 

m([-x,T-x\), ifT-W<x<0, 

m([-x,W]), if -W<x<T-W, 

Var[5(x)}, 

T 

fi(x)v(dx), 



w 

T 



f 

J-W 



X 



X(x)(u) 



f 2 (x)v(dx), 

- expectation of the distribution of claim sizes, if it exists, 

- variance of the distribution of claim sizes, if it exists, 
: D([-W,T}) -»• Rl '^ defined as 

- x(T — u) — x((—u)—), 



-J 

J[0,W] 



E[ X {N™){u)]m{du), 
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d 2 = [ [ Cov[x(N QO )(u),x{N 00 )(v)}m(du)rh(dv). 
J[0,W] J{0,W] 
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